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Dilute frozen solutions of the single molecule magnet Ni4 (5* = 4) have been studied using high 
frequency D-band (130 GHz) EPR. Despite the random orientation of the molecules, well defined 
EPR absorption peaks are observed, due to the strong variation of the splittings between the dif- 
ferent spin-states on magnetic field. Temperature dependent studies above 4 K and comparison 
with simulations enable identification of the spin transitions and determination of the Hamiltonian 
parameters. The latter are found to be close to those of Ni4 single crystals. No echo was detected 
from Ni4 in pulsed experiments, which sets an upper bound of about 50 ns on the spin coherence 
time. 



Single molecules magnets (SMMs) have been suggested 
as candidates for qubits in quantum processors^. How- 
ever, the coherence time (T2) of high-spin molecules has 
not yet been determined. Quantum tunneling of magne- 
tization (QTM) has been widely studied in SMMs^i^ii. 
A recent focus is on coherent QTM in which the tunnel- 
ing rates may be faster than the rate of decoherence^i^ii. 
Due to inhomogeneous broadening, only a lower bound 
(« 0.5 ns on most of the measured SMMs) of T2 can 
be extracted. Furthermore, in SMM crystals the strong 
dipolar interactions between the molecules (separated by 
only 1 nm) is expected to drastically reduce the coher- 
ence time. In order to determine the latter [e.g. through 
spin echo) and to coherently manipulate the magneti- 
zation with microwave pulses in SMMs, it thus may be 
necessary to work with dilute ensembles of molecules. 

SMMs may be diluted in solvents. The drawback of 
such a method is that the molecules are randomly ori- 
ented. However, magnetometry^ and EPR experiments*^ 
on randomly oriented powder have been succesfuUy per- 
formed. Only recently have dilute frozen solutions been 
studied by high frequency (HF) techniques. Frequency 
domain spectroscopy showed that the zero field splitting 
(ZFS) parameters of Mni2 clusters diluted in a frozen 
solvent were almost the same as in the crystal formic. 
Even more interesting, standard microwave pulse exper- 
iments have determined a long intrinsic coherence time 
(few microseconds at T = 2 K) in doped antiferromag- 
netic wheels (low-spin)ii. 

In this paper we present continuous wave (cw) HF EPR 
experiments performed on dilute frozen solutions of the 
SMM Ni4 {S = 4) as a function of temperature between 
6 and 30 K. Well defined transitions are observed in the 
0-7 Tesla range at / = 130 GHz. Simulations using 
ZFS parameters close to the ones determined from single 
crystalsi^ and D and g strains qualitatively reproduce 
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Figure 1: (Color online) EPR spectrum of a 10 mg/mL Ni4 
solution at / = 130 GHz and T = 5.5 K. The dashed and solid 
lines are respectively the absorption and dispersion channels 
(a 12 Gpp field modulation is used). The spurious signal at 
low field was also observed without Ni4 molecules in the 1:1 
mixture of toluene and dichloromethane. 

the spectra measured. No spin echo from the Ni4 clusters 
was observed in pulse experiments, which sets an upper 
limit for T2 of about 50 ns. 

[Ni(hmp)(dmb)Cl]4, henceforth referred to as Ni4, is a 
particularly clean SMM with no solvate molecules present 
in its crystal phase and only 1% (natural abundance) 
of nuclear spins on the transition metal sitesi^i^^. This 
results in narrower EPR peaks than in many SMMsi^. 
The spin Hamiltonian of Ni4 is to first approximation: 

n = DSl + BSi + C{SX + St) + iLsH.g.S, (1) 

where the first term is the uniaxial anisotropy, the second 
and third terms are 4*''-order anisotropy terms, and the 
last term is the Zeeman energy. The 5 = 4 ground state 
of the molecule at low temperature is the consequence 
of ferromagnetic exchange interactions between the four 
Ni^^ (5 = 1) ions. The uniaxial anisotropy leads to a large 
energy barrier « 12 K to magnetization reversal 
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between states of projection Sz = ±4 along the easy axis 
of the molecule. The ZFS parameters have been mea- 
sured on single crystals using HF EPRi^. In particular, 
the strong C-term (w 10 MHz) explains the fast tunnel- 
ing observed at zero fieldi^, and a recent analysis shows 
that the fourth order terms in Eq. [T] may result from a 
finite ratio of the (second order) single ion anisotropics 
to exchange constanti^. 

To prepare dilute solutions of Ni4, dried crys- 
tals were dissolved in a 1:1 mixture of toluene and 
dichloromethane. This glass was chosen due to sol- 
ubility restraints. ^HNMR spectra were collected 
for [Ni(hmp)(dmb)Cl]4 and the analogous complexes 
[Ni(hmp)(CH30H)Cl]4 and [Ni(hmp)(CD30D)Cl]4, 
as well as the free ligand 2-hydroxymethylpyridinc 
(hmpH). The proton signals for the bound hmp^ 
ligand were positively identified through compari- 
son of the [Ni(hmp)(CH3 0H)Cl]4 and deuterated 
[Ni(hmp)(CD30D)Cl]4 complexes. Comparison of 
the [Ni(hmp)(dmb)Cl]4 and free hmpH spectra con- 
firmed that no free ligand is present, and thus, the 
[Ni(hnip)(dmb)Cl]4 complex is intact in solution. Most 
of our EPR experiments were carried out on a 10 mg/mL 
solution, corresponding to an average distance of 6 nm 
between molecules. 

D-band (130 GHz) EPR spectra were acquired on a 
spectrometer capable of working in cw and pulse modes. 
The quadrature detection microwave bridge was designed 
and built by HF EPR Instruments, Inc. (V. Krymov, 
New York). During the acquisition in cw mode, a field 
modulation of 12 Gpp is generated at 100 kHz by a coil 
surrounding the cavity, and the dispersion and absorp- 
tion channels^i are recorded as a function of the exter- 
nal field, generated by a 7 T superconducting magnet. 
Fig. [1] shows the result of such an experiment performed 
at 5.5 K on a 10 mg/mL Ni4 solution contained in a 
quartz capillary tube (active volume 0.2 /iL). Large ab- 
sorption peaks (labeled a to C, with decreasing field) can 
be observed above 2.5 T. They are attributed to transi- 
tions between spin-states of the dissolved Ni4 molecules. 
The narrow signal around 4.6 T {g k, 2) and the spurious 
absorption at low field were also measured without Ni4 
clusters in the solvent, in contrast to the high field peaks. 
These do not affect the high field part of the spectrum, 
and we will not discuss them further. 

To address the origin of the transitions observed in 
Fig. [U we discuss the dependence of the EPR signal on 
temperature and the comparison with simulations, shown 
in Fig. [2^ and b respectively. The latter were performed 
using the EasySpin toolbosi^, which allows calculation 
of EPR spectraF^ and provides the visualisation tools 
used in Fig. [3l 2°x2° cells in one octant of the upper 
hemisphere were used to perform the powder average. 
First, we used the ZFS parameters of Ni4 single crystal 
(D = -0.748(5) K,B = -6x10"^ K, C = ±2.9x 10"^ K, 
9x = 9y = 2.23, gz = 2.3) ^^i^^ , and a convolution by a 
gaussian of width 0.1 T (dashed curve on Fig.l^b at 6 K). 
The position of the high field peaks arc slightly offset with 
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Figure 2: (Color online) a) EPR signal versus temperature, 
b) The dashed curve is the simulation at 6 K using the ZFS 
parameters of the single crystal (see text) and a 0.1 T gaus- 
sian broadening. The continuous curves have been calculated 
using the parameters indicated in the inset, that include D 
and g strains. 



respect to the data, and a large absorption peak close to 
zero field is calculated, but not detected. 

To improve the position of the calculated peaks, 
we used slightly different D and B parameters {D = 
-0.746(5) K, B = -9 X 10"^ K), and introduced D and 
g strains [au =Q-ll\D\ and ag = 0.02g), which wash out 
the low field transitions and lead to an overall better fit 
to the spectrum. It is known that D and g strains influ- 
ence the lincshapcs in SMM single crystals^. Random 
dipolar fields, also expected to play a role, should be less 
important in dilute solutions due to the larger average 
distance between neighboring molecules. The simulations 
reproduce rather well the dependence of the experimen- 
tal data on temperature. In particular, it can be seen 
that, on top of the global decrease of the absorption, the 
simulations reproduce the extinction of the peaks a and 
/3 as the temperature increases. This is expected from 
Boltzmann statistics, since the difference in population 
between the lowest energy levels decreases exponentially 
with temperature (note that / = 130 GHz corresponds to 
T ~ 6.2 K). Therefore, a is associated with a transition 
from the ground state to the first excited state. How- 
ever, we were unable to reproduce the lineshapes and the 
very intense e peak. We explored the possibility of hav- 
ing two Ni4 species (dissolved and prccipitatcdiSi) in the 
solution without success. We believe that our dilute so- 
lutions are free of precipitate, and it was checked after 



3 




Figure 3: (Color online) a) Energy levels for 3 different angles 
between the applied field and the molecular easy axis. The 
transitions corresponding to f — 130 GHz are shown. The 
stronger the red color, the larger the matrix element of the 
transition, b) Orientations contributing to transitions a and 
e {cf. Fig. [T]and[2|. Note that in a) and b), all allowed 
transitions are displayed without regard to the temperature. 

each experimental run that the dilute solution contained 
in the quartz tube was stiU clear. Increase of the D and 
g strains does not improve the fit either. The values used 
in the simulations are already larger than the ones fitted 
on single crystals^S by a factor 2 to 10. This probably 
results from the variation in the molecular environment 
of the SMM clusters in the frozen solution. It is known, 
6.17. in Mni2 crystals, that this environment plays a crit- 
ical role^i. In solution, the strain could also depart from 
axially symmetric. i?-strain can not be excluded, as sim- 
ulations performed with cte = 0.05 K slightly broaden 
peaks a to 7. 

To better understand the nature of the transitions in- 
volved in the EPR spectrum. Fig. [3^ shows the energy 
levels as a function of the angle 9 between the easy axis 
of the molecules and the magnetic field (the angle in 
the hard plane is set to 0°). The well-known diagram 
is obtained for 6 = 0°, with the low lying transitions 
at low field. In our experimental case, the ground state 
transition (m = —4 — > —3) at low field is not observed 
because the ZFS of Ni4 is about 140 GHz, slightly larger 
than our working frequency. When 6 = 90°, the ground 
state transition is on the contrary at the highest field {a 
in Fig. [1] and [2). For an intermediate angle, the tran- 
sitions can also be determined by a full diagonalization 
of the Hamiltonian (Eq. [1]). As a result of the rapid 
variation of the splittings with the amplitude and the 
direction of the applied field, certain subpopulations of 



molecules are selected, and well defined absorption peaks 
are observed. Fig. [3)d, which shows the orientations con- 
tributing to peaks a and e, validates this picture. The 
influence of the C-term can be seen for the peak a (vari- 
ation between = 0° and </> = 45°), for which the easy 
axis of the contributing molecules almost lie in the plane 
= 90° (the band near = 0° is a highly excited transi- 
tion which does not contribute at low temperature). 

Finally, echo experiments (7r/2 — r — tt — t— echo high 
power pulse sequences) were conducted at 5.5 K on sev- 
eral dilutions of Ni4 (down to 0.5 mg/mL) and on the 
pure solvent. We did not observe any echo from the Ni4 
molecules. Let us briefly discuss this result. First, the 
sensitivity of the experimental setup is sufficient to detect 
echoes with T2 > 50 ns in 200 G broad lines from micro- 
molar solutions (i.e. 10^"'^ spins) of 5 = 1/2 compounds. 
Despite the random orientation of the Ni4 molecules and 
the broad absorption lines (w2 to 4 kG), the cw EPR 
data shows that a small («5%) but non negligible frac- 
tion of the molecules contribute to the signal. Taking 
into account these factors, it is found that there are still 
about 30 times more spins than needed to detect an echo 
if T2 > 50 ns for the 10 mg/mL Ni4 solution. Second, 
the random orientation makes the intermolecular interac- 
tion weak because the effective dilution of the molecules 
adressed by the microwaves is large. Third, the coupling 
to nuclear spins can yield coherence times larger than 
50 naii. It is not clear whether the coherence time of 
spin-states oscillations between excited levels (as /? to 
could be much shorter than the one between the ground 
state and the first excited states (as a). However, an 
estimation of the spin-lattice relaxation time from direct 
spin-phonon process^^ using a reasonable value for the 
transverse speed of sound in the frozen solvent (1000 m/s) 
yields Ti = 40 ns at / = 130 GHz and 5.5 K (it increases 
by a factor 2 at K). This upper bound for T2 could 
explain the absence of echo from Ni4 in our experiments. 
It is interesting to note that Ti is a strong function of 
frequency (~ 1//^)- For instance, at 10 GHz and 1.5 K 
we estimate Ti = 60 /xs. 

In conclusion, we have studied dilute frozen solutions of 
the SMM Ni4 by HF EPR. The clusters have been found 
to be stable in solution. Despite the random orientation 
of the anisotropy axes, it is possible to measure EPR 
spectra with well defined transitions and to deduce that 
the ZFS parameters of the dissolved molecules are close 
to the ones in the crystal phase. Thus, the behavior of 
an isolated Ni4 cluster is essentially the same as that of 
a single crystal. 

This work was suppported by NSF-NIRT Grant No. 
DMR-0506946 (ADK) and by NIH Grant No. GM075920 
(GJG). 
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